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We analyze operational risk in terms of a spin glass model. Several regimes are investigated, as 
a functions of the parameters that characterize the dynamics. The system is found to be robust 
against variations of these parameters. We unveil the presence of limit cycles and scrutinize the 
features of the asymptotic state. 
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I. INTRODUCTION 

The physics of complex systems has proven to be a fer- 
tile arena of investigation in very diverse areas of natural 
sciences. The versatility of the techniques of analysis and 
the wide applicability of the main underlying ideas make 
this subject interdisciplinary: there are applications in 
physics, mathematics, statistics, cybernetics, chemistry, 
biology, as well as economics and social sciences, this list 
being far from complete. 

The guiding principles in the study of complex systems 
are philosophically intriguing, and lie in the idea that 
one can learn simple lessons from complexity [1]. This 
requires the definition of an appropriate level of descrip- 
tion at which the main features of the system investigated 
can be described in terms of a given number of suitable 
variables. The very definition of these variables and their 
interactions design a convenient "minimal model," whose 
dynamics, equilibrium features and correlations can be 
investigated. 

Spin glasses play an important role in this context 
[2, 3]. A spin glass can be thought of as an ensemble of 
interacting spins, whose bonds are randomly distributed. 
The magnetic ordering that emerges resembles the posi- 
tional ordering of a physical glass and is characterized by 
frustrated interactions. Like all complex systems, a spin 
glass is composed of interconnected parts (spins) and as a 
whole exhibits properties that arc far from being obvious 
if one looks at the properties of the individual parts. 

In this article we shall apply a spin glass model to 
study a financial system [4, 5]. We shall seek an appro- 
priate level of description at which capital losses can be 
schematized as interacting dichotomic variables (spins). 
We shall focus in particular on "Operational Risk" (OR) 
[6, 7], defined by the The New Basel Capital Accord 
(Basel II) [8] as "the risk of [money] loss resulting from 
inadequate or failed internal processes, people and sys- 
tems or from external events". OR management [9] is 
a subcategory of risk management, that branch of econ- 
omy whose purpose is the identification of the possible 
sources of money losses and eventually the definition of a 
strategy for avoiding them. According to the regulations 
set by the "Basic Indicator Approach" proposed in Basel 



II, banks have to set aside 15% of their total capital, in 

order to face operational losses: since this is a quenched 
capital, banks share a keen interest in order to keep this 
percentage as small as possible. 

Basel II also proposes a slightly refined approach (the 
"Standardized Approach") that defines eight different 
"business lines" [16] on which losses can occur, by as- 
sociating to each of them a constant risk coefficient that 
can vary between 12% and 18%, so that the total amount 
of money to be stored by a bank is obtained by a weighted 
sum. The Basel accord gives to each bank the possibility 
to propose an internally developed "Advanced Measure- 
ment Approach" to estimate the capital to set aside. Let 
us point out that, while the Basic Indicator Approach 
and the Standardized Approach arc aimed only to the 
estimation of this capital, the Advanced Measurement 
Approaches are particularly attractive for a bank because 
they may be oriented to the management of OR. Every 
Advanced Measurement Approach is required to classify 
the operational losses in the 8 afore-mentioned business 
lines and in 7 "event types" [17]. There are therefore 
56 combinations (8 business lines x 7 event types) to be 
implemented. The objective of our study is i) to accu- 
rately characterize these different losses, in order to ii) 
understand the main features of their dynamics and in- 
teractions, and if possible iii) minimize them. 

In order to analyze the distribution of these losses, we 
shall model them via N{= 56) interacting variables, as 
proposed by Khiin and Neu [10] and elaborated by several 
authors in Refs. [11-13]. This sets the afore-mentioned 
level of description at which techniques of spin glasses 
and complex systems can be applied. An efficient perfor- 
mance of this model requires knowledge of the frequency 
and value of each loss in the N different channels and a 
distribution of such events. Because Basel II is a rather 
recent treaty, real data are often insufficient. For exam- 
ple, some of the 8x7 channels are often completely empty, 
so that rare events cannot be included. Also, some busi- 
ness lines (e.g. Retail Banking) are often very populated, 
so that they tend to "overwhelm" other lines when a sta- 
tistical analysis is performed. Finally, the populations 
can be very different for different banks, and can show 
significant variations according to, e.g., bank size and ge- 
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ographical location. To focus, say, only on large banks, 
would not be meaningful, because one could aim at offer- 
ing a service to many small banks (say, minimize below 
15% the percentage of the total capital to be put aside), 
by taking into account their peculiarities and characteris- 
tics in a given local context. To add to these difficulties, 
banks are often reluctant to disclose the details of their 
losses, for obvious reasons. We shall therefore opt for a 
democratic approach, by setting the parameters of our 
model to reasonable values, valid for a rather wide range 
of banks, and checking that our results are stable against 
small variations of these parameters. 



II. THE MODEL 

The model we shall study consists of N dichotomic 
interacting variables that evolve in time. These variables 
schematize the 8x7 interacting processes defined in the 
Introduction. The model was introduced in Ref. [10], 
elaborated in [11] and [12], where processes was proven to 
be dependent on each other, and further investigated by 
Khiin and Anand in Ref. [13], where correlations among 
processes were further explored. We shall adopt in the 
following a slightly different notation from the preceding 
articles. See Eqs. (9) and (10) in the following. 

Let s{t) be the A^-dimensional vector of the dichotomic 
variables and the coupling constant between Sj and 
Sj. Time is discretized t e N and Si{t) G {— 1,+1}, 
with i € {1, • • • tN}. The value —1 denotes a nmning 
process (no losses), while the value -|-1 a broken down 
process, that does not run anymore. In such a way, we 
are completely neglecting the amount of the loss and are 
only recording whether a loss occurs or not. In particular, 
since banks record only losses exceeding a given threshold 
[18], we can interpret Si{t) = +1 as the occurrence at 
time t of a loss equal to or larger than this threshold in 
the i-th channel. 

We assume that the time evolution is governed by the 
equation 



Si{t + 1) = sign 



N 



(1) 



where 9i are characteristic constants, which we call sup- 
ports, and ^i{t) are Gaussian noises 

(ei(f))=0 and {^i{t)^^{s)) = 5i,j5t,sTl (2) 

Here the brackets (•) denote the average over the noise 
and 8i^j is Kronecker's delta. 

In order to understand the role of the parameters that 
appear in the evolution equation, let us decouple the pro- 
cesses, by setting Jjj = V so that Eq. (1) reads 



The probability that Si(t + 1) = +1 is easily calculated 

Am 



Pi = P{si(t+1) = + 1} 

that, after integration, yields 
1 



e (4) 



Pi 



1 - erf 



(5) 



where erf(2:) = 27r~^/^ J^^ e^*^/^ dt is the error function. 
(Notice that Si{t + 1) is decoupled even from Si{t) in this 
simple example.) We see that a given process Sj runs 
after the time increment t — >■ i -I- 1 if it is "supported" 
by a suitable environment, that yields an average sup- 
port 6i{> 0). This elucidates the role of the supports for 
noninteracting variables: when nonvanishing couplings 
appear, the evolution is given by (1), so that a process 
keeps running over the time increment t — ^ t -|- 1 if it 
is suitably supported by the argument of the sign func- 
tion, which is a combination of its average support and 
the action of the other processes. Observe that the sign 
of the couplings Jij is not defined and that the average 
support 0i is presumably defined by a process-dependent 
environment: a positive value 6*, > denotes a certain 
"care" of the bank that process i keeps running. 

The evolution equation can be given physical meaning 
by defining the following function of the configurations: 

^ N N 

^ (*) = - 9 E (*) (*) + E ^'^^ (6) 



i=l 



If the couplings are randomly distributed, the above 
"Hamiltonian" describes a spin glass in a site-dependent 
magnetic field. If at each time step only one variable 
evolves (the so called asynchronous updating), one can 
prove after some algebra that the function AH{t) = 
H{t + 1) — H{t) is negative only if the following condition 
holds: 



sign 



sign 



E'Jij ~\~ Jji 



At) 



(7) 

This means that, under the previous condition, H{t) is 
a monotonously decreasing function of the configurations 
as the system evolves according to Eq. (1). This equation 
of motion can be independently derived from Eq. (6): 



Sk (i + 1) = sign 



dH 
dsk 



it) 



(8) 



s^{t + l)=sign i-e^ + m)]- 



(3) 



If H is viewed as a potential, Sk aligns to the force derived 
from this potential. In the case of synchronous updating, 
which is our case, one numerically observes that AH{t) 
is negative if the coupling constants are averaged over. 
Observe that the condition (7) is not entirely obvious in 
our context, where, say, a very lossy and highly populated 
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channel might influence more a rarely populated channel 

than vice versa. 

In order to connect our results to those in [13], notice 
that if we bijectively map the variables s £ {— 1,+1} into 
the variables t] G {0, +1}: 



s,it) = 277,(t)-le {-!,+!}, 
Vi{t) = I [Si{t) + 1] G {0, 1}, 



(9) 



Eq. (1) reads 



Vi{t + 1) = -sign 



= e 



= e 



i 

j 



+ 



(10) 



where Q is the Heaviside step function, and Jij = 2Jij 
and 9i=ei + Y,j Jij- 

III. ROBUSTNESS OF THE MODEL 

Before we start investigating the model introduced in 
the preceding section, it is useful to study its robustness 
when the parameters that characterize it are changed. 
This will also help us get a feeling of the way the dynam- 
ics is implemented and the different roles of the parame- 
ters. 

In order to study the evolution over T time steps, let us 
define the vector of integers z{T) = {z\(T), . . . , zn{T)), 
with 



1), 



(11) 



that records the positive (+1) outcomes of the s = 
(si, . . . , spf) variables during the whole evolution T. Note 
that the integers Zi(T) are bounded by < Zi(T) < T. 
We ask the following question. Let G N-*^ be a given 
vector (for instance the data set obtained by accumulat- 
ing the losses of a given bank during a certain number 
of years). How should we choose the parameters of the 
model dynamics in such a way that the evolution of the 
system yields at time T? The simplest way to do this 
[13] is to assume that all variables are independent and 
define the probability that a variable takes the value 1 
according to 



Pi 



(12) 



In this way, we are assuming that these probabilities are 
constant in time. We then set in Eq. (1) all couplings 



Jij = (independent variables), and from (5) the sup- 
ports 



^i = -y2^erfinv(2pi-l). 



(13) 



where erfinv is the inverse of the error function erf, and 

let the system evolve for T = time steps. In this 

way we recover the given z^ on average. Indeed, 



{zi{T)) = Tpi = 



(14) 



We now introduce interactions among processes and 
consider the coupling constants as Gaussian random vari- 
ables with average and width 



and 



(15) 



respectively, where the bar denotes the average over the 
distribution of the Jjj 's. This is the familiar quenched 
picture: couplings change over a slow time scale and can 
be considered constant during the whole evolution, up to 
time T. The average must be performed at time T over 
several evolutions, each with a different realization of the 
couplings. 



TABLE I: Data set of accumulated losses Zt and correspond- 
ing supports 9 ior N = 4. See upper panel of Fig. 1. They 
have mean fj,e = 0.7017 and standard deviation ae = 0.8372. 
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1.2688 


1.5287 



TABLE 11: Data set of accumulated losses z» and correspond- 
ing supports 6 for A'^ — 10. See lower panel of Fig. 1. They 
have mean ne = 0.9516 and standard deviation ae = 0.5249. 



4637 2994 
90 30 



1648 
90 



239 
60 



209 
3 



e 0.0456 0.2630 0.4874 0.9893 1.0174 
1.1833 1.3744 1.1833 1.2566 1.7162 



We start from the initial configuration z(0) = (all 
processes running) [19] and numerically evaluate the dis- 
tance 



mT))-z4 

Tx/N 



(16) 



where the averages are defined in Eqs. (2) and (15) and 
the norm is the Euclidean one, jlajl = iY^iCif)^^^ ■ (Since 
< {zi{T)), Zi:^i < T, the normalization is such that 
< ^ < 1.) We look at two cases: N = A, with the 
supports shown in Table I, and N = 10, with the sup- 
ports shown in Table II. The results are shown in Fig. 
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In order to test the robustness of this result, let us add 
a constant external field a = {a, . . . ,a) to the supports 



(17) 



We shall consider a small (a — 0.2) and a large {a — 2) 
perturbation of the supports 6. The presence of a intro- 
duces a bias (a systematic error) in the dynamics, and 
we ask whether it is possible to compensate for this er- 
ror by suitably tuning the interactions. The dashed line 
[a — 0.2) and the thick solid line (a — 2) in Fig. 1 an- 
swer this question: the distance is large in absence of 
spin interactions (tJj = 0). However, for larger ctj, the 
distance reaches a minimum because some (random) cou- 
plings correct the evolution of the system in the desired 
direction. 

These results can be read in different ways. Interac- 
tions among processes can partially compensate for an 
erroneous choice of the supports. The very presence of 
a minimum (characterizing "optimal" couplings) is an 
interesting feature of the model. Vice versa, strong in- 
teractions among processes can lead the dynamics astray, 
nullifying the stabilizing (beneficent) action of the sup- 
ports. The interplay among the different parameters will 
contribute to make the (highly nonlinear) dynamics a 
very interesting one. 



FIG. 1: (Color online) Distance 5 in Eq. (16) between the 
reference and final configurations vs standard deviation aj 
of the distribution of the coupling Jij. Upper panel A'^ = 4; 
lower panel A'' — 10. Thin solid line: vanishing external field 
a. The minimum is at aj — 0, when all couplings vanish, and 
the final configuration is identical to the initial one; as ctj 
increases, the distance monotonically increases. Dashed fine: 
a = 0.2; thick solid line a = 2. In both cases the distance 
is no longer zero at crj = because of the bias introduced 
by the field. As aj increases, the spin interaction partially 
compensate for the external bias, so that the distance reaches 
a minimum and then increases. 



1 (thin solid line a = 0). The evolution yields a final 
vector z{T) very close to z* for aj ~ 0, as expected. As 
CTj increases, the spins interact, causing the distance S to 
increase monotonically. The result is qualitatively simi- 
lar for both values of N. We observe that the evolution 
yields a vector that is rather close to i** for sufficiently 
small couplings. 



IV. ANALYTICAL APPROXIMATION 

We are now ready to tackle the problem in full general- 
ity. Let us start by considering a one-step evolution, that 
can be analytically averaged. It is convenient to rewrite 
Eq. (1) in the form 



Si(t+l) 26 



^J,,s,(^)-0, + e^(^) 



- 1- (18) 



We now substitute the Heaviside function as the integral 
of its Laplace transform (l/z) 



I f dz 



along a vertical line B just at the right of the imaginary 
axis (Bromwich's path). 

We take the average over the noise and the couplings 



n 



dJi 



^2^ 



=e « 



27rcr| 
I 



(20) 



Gaussian integrations yield 
1 



On the Bromwich path, the integration variable takes 



(21) 



5 



the form z = 7 + iy with 7 > 0; the exphcit form of the 
previous integral is then 



X [cos [jyk - ydi] + isin [^yk - y6i\\ - 1, 



where 



(22) 
(23) 



In the hmit 7 — )• 0, the first contribution yields 
1 



Ay 



7r J-00 y 



e cos [7?;fc - y^i] 



(24) 



in which we have used the well-known formula for distri- 
butions 



1 



X ~ {xq ± ie) X — Xq 



±iTrS(x-XQ), (25) 



with P denoting the Cauchy principal value. Due to 
symmetry considerations, the previous integral yields a 
contribution 1/2 to the total integral. The remaining 
part does not have a pole, so we can set 7 = in the 
integral, obtaining 



if* 

J -00 

Summing up, 



sin(-y6li) 



e 2« "5 



-erf 



(s.(t + l)) 



erf 



2{al 



(26) 



(27) 



It is clear that by iterating this expression over the whole 
evolution [0, T] one makes a mistake, because the average 
over the couplings is taken at every time step. 

Another alternative would be to compute the noise- 
averaged single-step evolution 



_ 5f 

e ^ 



27rcr| 



xsign 



N 



JijSj(t) -9i+ 



erf 



(28) 



and average over the realizations of the couplings at the 
final step. This also entails an error, because we are 
forcing a Markovian evolution. 




FIG. 2: (Color online) Evolution of the magnetization. Com- 
parison between the numerical evolution (1), and the approx- 
imations (27) and (28). Notice that m(0) = 0.5. 




FIG. 3: (Color online) Evolution of the quantity z{t) = 
'^iZi{t), see Eq. (11). Comparison between the numerical 
evolution (1), and the approximations (27) and (28). 



We shall now compare the expressions (27) and (28) 
to the numerical evolution (1). In particular we will look 
at the quantity 



(29) 



which is simply the sum of the mean components of the 
spin vector s and has the (appealing) physical meaning of 
average magnetization. This quantity will be thoroughly 
analyzed in the remaining part of this article. 

The three evolutions in Fig. 2 are similar, in that the 
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magnetization tends in all cases to an asymptotic value 
m — lim m{t). (30) 

t— f oo 

These asymptotic values coincide within 10%. Notice, 
however, that this error (linearly) accumulates in Eq. 
(11), as shown in Fig. 3. Moreover, both approxima- 
tions (27) and (28) are unable to describe the features of 
the short-time behavior. For this reason, we shall numer- 
ically evaluate the evolution of the magnetization (29), 
by using (1). 

V. ASYMPTOTIC MAGNETIZATION 

The magnetization (29) and its asymptotic value (30) 
clearly depend on the details of the model. However, 
the initial spin configuration has practically no influence, 
as shown in Fig. 4, that displays the evolution of m(t) 
starting from three different initial configurations: one in 
which all processes are initially running, Si(0) = —1 Vi, 
one in which they are all broken down, Si{0) — 1 Vi, and 
one in which they are randomly distributed. We took 
N = 10, = aj = 1 and the supports (randomly) dis- 
tributed in [-1, 1], namely = (0.3115, -0.9286, 0.6983, 
0.8680, 0.3575, 0.5155, 0.4863, -0.2155, 0.3110, -0.6576), 
with mean 

1 ^ 

^« = ]^E^^ (31) 

1=1 

equal to 0.1746 and standard deviation ag = 0.5870. 
Clearly, the system quickly forgets its starting point, and 
the magnetization always converges to a common state 
that does not depends on the initial configuration. For 
this reason in the following analysis we will always con- 
sider an initial configuration in which all spins are down. 

By contrast, the magnetization strongly depends on 
the supports and on the standard deviations of the dis- 
tributions of the couplings and the noise. When and 
aj are very large, the asymptotic magnetization always 
tends to 0, independently of the supports. This can be 
seen in all the figures that follow and is in agreement 
with intuition: when noise is large and the distribution 
of couplings has a large standard deviation, the role of 
supports becomes negligible in Eq. (1), and the spins take 
the values ±1 with equal probability. 

Let us first consider the case of = 5 spins and 
large positive supports in the interval [4, 6] , namely 9 = 
(4.1951,4.5570,5.0938,5.9150,5.9298) (fie = 5.1381 and 
ag = 0.7841). We first consider the non-interacting case 
(a J = 0). As CTj increases, the spins can make transi- 
tions to the -fl state and the asymptotic magnetization 
increases. Since all supports are positive, the standard 
deviation of the noise must exceed a certain thresh- 
old before transitions can take place. See solid line in 
Fig. 5. On the other hand, for strongly interacting spins 
a J = 1.6, the effect of noise is essentially negligible, see 
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FIG. 4: (Color online) Time evolution of the magnetization of 
the system for three different initial configurations: all spins 
up (dots), all spins down (stars) and randomly oriented spins 
(squares). The couplings are quenched random variables, cho- 
sen according to (15), while the supports 9i £ [—1, 1] with 
fio = 0.1746 and ag = 0.5870. 
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FIG. 5: (Color online) Asymptotic magnetization vs a^. The 
couplings are quenched random variables, chosen according 
to (15), while the supports 8i £ [4,6] with jj,e — 5.1381 and 
ag — 0.7841. Solid line: aj = (no couplings); dashed line: 
aj = 1.6. 



dashed line in Fig. 5. Clearly, both lines tend to as 
CTj — >■ oo (not shown). 

We now perform the same analysis by exchanging the 
role of the two standard deviations. Figure 6 displays the 
behavior of the same system with large positive supports 
in [4, 6], as aj varies. The solid line is the evolution in the 
noiseless case, while the dotted line is the evolution for 
CTj = 1.6. The behavior is qualitatively similar to that in 
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FIG. 6: (Color online) Asymptotic magnetization vs crj. The 
couplings are random quenched variables chosen according 
to (15) while the supports 6i £ [4,6] with /ig = 5.1381 and 
ae — 0.7841. Solid line: = (noiseless case); dashed line 
0-5 = 1.6. 




FIG. 8: (Color online) Asymptotic magnetization vs aj and 
The couplings are quenched random variables, chosen 
according to (15), while the supports 6i G [—1, 1] with = 
0.3484 and ag — 0.6974. The magnetization tends to zero as 
and/or aj become large. 
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FIG. 7: (Color online) Asymptotic magnetization vs aj and 
(75. The couplings are quenched random variables, chosen 
according to (15), while the supports 6i £ [4,6] with ne = 
5.1381 and erg = 0.7841. The magnetization tends to zero as 
(jj and/or aj become large. 



Fig. 5, but quantitatively different. We shall come back 
to this point in the next section. 

These conclusions are summarized in Fig. 7. Notice the 
asymmetry O aj. Observe that the magnetization 
tends to vanish as cr^ and/or aj become large. Notice also 
that due to the symmetry of the problem, the situation 
with large negative supports would be identical (modulo 
a reflection m o —m). 

The scenario changes when the supports are ran- 



domly chosen in the interval [—1,1], e.g. = 
(-0.6848,0.9412,0.9143,-0.0292,0.6006) {^le = 0.3484 
and <jg — 0.6974). Figure 8 shows the asymptotic mag- 
netization when (7j and cr^ vary. As one can see, the 
magnetization (quickly) tends to a minimum and then 
slowly tends to its asymptotic value 0. 

The plot in Fig. 9 illustrates this aspect in more de- 
tail. Let us consider the case with no couplings, aj = 0, 
and look at the asymptotic mean value of two 'particu- 
lar spins: one with a large (positive) support Q = 0.9412 
and one with a small (negative) support d = —0.0292. 
When 1^1 is small, a small noise is sufficient to drive 
the mean value of the spin to zero. On the other hand, 
a spin with a large support is rather insensitive to small 
values of noise cr^, which should therefore be sufRciently 
high in order to drive it towards 0. By summing these 
different contributions to the asymptotic magnetization 
one obtains the minimum in Fig. 8. More in general, 
the behavior of the magnetization for small values of 
and/or cj strongly depends on the value of 6'i's, while 
(obviously) as the standard deviations increase, all vari- 
ables tend to the same asymptotic value, and this yields 
the flat part of the graph. This qualitatively explains all 
the features of Fig. 8. We remind here that Si — +1 rep- 
resents a loss in channel i, while = — 1 means that no 
loss occurs. 



VI. ASYMPTOTIC BEHAVIOUR 

In this section, at variance with the previous ones, we 
will consider couplings with a nonvanishing mean. How- 
ever, we will appropriately scale the mean and standard 
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FIG. 9: (Color online) Asymptotic average value of two par- 
ticular spins vs (Tj. The two spins have very different supports 
9. The couplings are quenched random variables, chosen ac- 
cording to (15). When \6\ is small the convergence is very 
fast, not so when it is larger. By summing curves like those 
in this figure one obtains the minimum that appears in Fig. 8. 



deviation of the couplings in term of the lattice size N 
and the mean support fj,g in Eq. (31) in order to disclose 
an interesting behavior. We take 



P(C.(t))=AA(0,a5), 



(32) 



where Af{p, cr) denotes the Gaussian probability density 
with mean /z and standard deviation a. The behavior 
of the asymptotic magnetization in the scaled variables 
is shown in Fig. 10 for N = 10. We take randomly dis- 
tributed Oi G [1,2], namely = (0.6294, 0.8116, -0.7460, 
0.8268, 0.2647, -0.8049, -0.4430, 0.0938, 0.9150, 0.9298) 
{ijg = 0.2477 and ag = 0.6917). Let us expHcitly observe 
that such value of introduces in (32) an antiferro- 
magnetic bias in the couplings. Remarkably, the graph 
displays a symmetry aj ^ as a consequence of the 
performed scaling. 

We now scrutinize an interesting feature of the evolu- 
tion equation (1). In the noiseless case this yields 



1 if J:f=lJ^JSJit)-0,>O, 



S^{t+^) 



-1 if Ef=i^uS, W-^^ <o. 



(33) 



If Jij = J < for alH, j, the previous equation reads 



S^(^+ 1) = 



+1 if m(t) < ^, 
-1 if m{t) > 



(34) 



Since TO(t) G {-1, -1 + ^, . . . , +1 - ^, +1}, all the sites 
such that -j^ < —1 {-jj^ > 1) are frozen in the —1 (-f-l) 



FIG. 10: (Color online) Asymptotic magnetization vs aj and 
CTj. The couplings are chosen according to (32), while the 
supports 9i G [1, 2] with fie = 0.2477 and ae = 0.6917. The 
magnetization tends to zero as and/or aj become large. 




FIG. 11: (Color online) Case in which all the supports satisfy 
the condition 9j/JN < Si, except the i-th one, that satisfies 
the condition Si < 9i/JN < ^2. The i-th spin flips, and the 
magnetization indefinitely oscillates between the values — 1 
and -1 + 2/7V. 



state. Let us suppose that this is valid for all spins ex- 
cept the i-th one, that can therefore flip. Only two pos- 
sible magnetizations are possible in this case: one with 

51 — +1, which we will call iS'2, and one with Si — — 1, 
which we will call Si < S2; if Si < 9i/JN < S2, Si oscil- 
lates between the up and down states. A very simple case 
is shown in Fig. 11: all the supports satisfy the condition 
dj/JN < Si except the i-th one, that satisfies the con- 
dition 5*1 < Oi/JN < S2- This means that 5i = — 1 and 

52 — —1 + 2/N. When the magnetization of the system 
is —1, the z-th spin flips, so the magnetization becomes 
— 1 + 2/N, then the spin flips again bringing the magne- 
tization back to the previous value. This oscillation goes 
on indefinitely. 

If — 2 spins satisfy the condition 9j/JN < Si and 
2 of them satisfy the condition Si < 9i < S2 a similar 
behavior occurs, but S2 — Si = A/N. Of course these os- 
cillations take place only if the interactions are negative. 
Incidentally, this explains why no limit cycle is observed 
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in Fig. 10. 

Summarizing, the above discussion implies that when 
the above-mentioned conditions hold, two different 
asymptotic magnetizations appear, one for even and the 
other one for odd times. Both tend to vanish when the 
standard deviations <tj,cf^ — >■ oo. We have thus unveiled 
the appearance in the system of an attractive limit cycle 
of period 2 [14, 15]. As will be clear from the following 
figures, the size of the cycle depends on the supports and 
shrinks for large values of the noise and/or the couplings. 
This yields important information about the typical val- 
ues of CTj and a J such that the evolving system does not 
feel any difference between even and odd times. 

Figure 12 displays the behavior of the asymptotic mag- 
netization for even and odd times and for two different 
realizations of the supports. In the upper panel 9 = 
(0.8636, 1.1562, 0.5200, 0.1114, 0.1399, 0.8885, 1.2400, 
0.2851, 0.8793, -0.0849) {fie = 0.5999 and erg = 0.4684), 
while in the lower panel 9 = (1.1576, 1.9706, 1.9572, 
1.4854, 1.8003, 1.1419, 1.4218, 1.9157, 1.7922, 1.9595) 
{fig = 1.6602 and ag = 0.3308). The values are averaged 
over several realizations of the couplings, distributed ac- 
cording to (32) . The main difference is the behavior near 
the origin: clearly, larger supports entail more inertia 
in the spin dynamics and a slower convergence to the 
asymptotic value. 

In Fig. 14 all supports are equal {ag = 0). In the left 
panel 9i — 9 = 0.1, Vz, while in the right panel 9i = 9 = 
0.8, Vi. As is to be expected, the manifolds appear to be 
more regular, in particular for larger 9. 

In Fig. 13 we focus on the role of the standard de- 
viation of the distribution of the supports. Again, the 
supports are quenched random variables, drawn from a 
Gaussian distribution with fig — 0.5. In the upper panel 
ag = 0.5 and in the lower panel ag — 1. The manifolds 
are similar and differ only for their behavior near the ori- 
gin: as stressed before, in that region the behavior of the 
magnetization strongly depends on the supports ^^'s (see 
Fig. 9). The two manifolds merge more easily for larger 
ag. 

Interestingly, a curious feature emerges in the case of 
uniform supports: Figure 15 shows the value of the stan- 
dard deviation at which the distance between even and 
odd magnetizations (radius of the limit cycle) is 1/4. In 
the upper panel a^ — aj; in the central panel = 0; 
in the lower panel aj = 0. A linear fit yields excellent 
accord with the numerical data. We offer no explanation 
for this linear behavior. 



VII. CONCLUSIONS 

We analyzed loss risks of banks by borrowing tech- 
niques of classical statistical mechanics. Losses were di- 
vided into N categories, each of which corresponds to 
a channel, schematized as a dichotomic variable (spin). 
This is paramount to assuming that one is only inter- 
ested in the occurrence of a loss and not in its amount. 



s 



s 




FIG. 12: (Color online) Asymptotic magnetization vs aj and 
CTj. The couplings are chosen according to (32), while the 
supports are extracted from a Gaussian distribution. Upper 
panel: fie ~ 0.5999 and ae — 0.4684; lower panel: fie = 
1.6602 and ag = 0.3308. A limit cycle appears: the two 
leaves represent the asymptotic magnetization for odd and 
even times. 



The basic idea is that a loss occurring in a channel can 
stimulate (inhibit) losses in other channels: this is mathe- 
matically described by assuming that loss events interact 
with each other by means of positive or negative cou- 
plings constants [10-13]. We assumed that the interac- 
tions are random, as in a glassy system. 

The time evolution of this model is governed by Eq. 
(1) and we studied the asymptotic behavior of the sys- 
tem for interesting values of the parameters. Our study 
brought to light some interesting features of the model. 
First of all, the system has a short memory: it quickly for- 
gets its initial spin configurations and the magnetization 
quickly converges towards its asymptotic value, which is 
determined by the supports and the parameters of the 
distribution function of noise and disorder. One should 
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FIG. 13: (Color online) Asymptotic magnetization vs aj and 
(75. The couplings are chosen according to (32), while the sup- 
ports are quenched random variables extracted from a Gaus- 
sian distribution with /le = 0.5. Upper panel ae = 0.5; bot- 
tom panel ae = 1. Like in Fig. 12, a limit cycle appears: the 
two leaves represent the asymptotic magnetization for odd 
and even times. 




FIG. 14: (Color online) Asymptotic magnetization vs aj and 
The couplings are chosen according to (32), while the 
supports are all equal: in the upper panel, 9i = 6 = 0.1, Vi, 
while in the lower panel 6i = 9 = 0.8, Vi. Again, a limit cycle 
appears, the two leaves representing the asymptotic magne- 
tization for odd and even times. The manifolds appear to be 
more regular than in Figs. 12 and 13. 



remember that the magnetization represents the average 
number of channels in which losses occur. 

Interesting features appear when the parameters are 
rescaled as in Eq. (32) and the supports are positive (we 
remind that a positive support means that a bank gives 
"care" for a source of potential losses in a given chan- 
nel) . For small values of the standard deviations of noise 
and disorder the system displays an oscillatory behav- 
ior (between even and odd time steps). This is a limit 
cycle [14, 15]. For larger values of these standard devi- 
ations, the system reaches a unique asymptotic value. 
This yields a scale on which the behavior of the sys- 
tem can be regarded as stable. In general, this scale 
depends on the parameters of the model. However, when 
the supports are identical, this scale linearly depends on 



the (common) support. We have offered no explanation 
for this curious numerical observation. 

In recent years, the scientific interest for economical 
issues has increased, thanks also to the contribution of 
advanced physical and mathematical techniques, such as 
stochastic processes, chaos theory and the theory of dis- 
ordered systems [4, 5]. This work is an effort towards 
the modelling of the operational risk of banks activities. 
Much additional work must be done, in particular in or- 
der to understand whether a glassy phase transition ap- 
pears. An additional point of interest, also in view of 
potential applications, would be to work towards the op- 
timization of some crucial parameters of the model, such 
as the supports. This would enable banks to understand 
in which direction their efforts should be directed in order 
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to minimize losses. 
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FIG. 15: (Color online) Values of the standard deviations at 
which the radius of the limit cycle (difference between magne- 
tization at odd and even times) is 1 /4. We look at three differ- 
ent directions: crj = (upper panel); cr^ = 0, (central panel); 
and (Tj = (lower panel). The linear fits <j{6) — a9 + b have 
slopes a = 0.45, 0.58, 0.82, respectively. Their correlation co- 
efficients p are bounded by > 0.9996 and the coefficient 6 
is compatible with 0, namely |6| < 8.6 x 10~*. 
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